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Abstract 

Recently Tsai et.al. (cond-mat/0406174) have used the renonnaliza- 
tion group approach to study strong couphng superconductors without 
assuming a broken symmetry phase. We use the hydrodynamic formu- 
lation to study the same problem with the same intention. We recast 
the electron-phonon + electron-electron problem in the hydrodynamic 
language and compute the one-particle electron Green function at finite 
temperature. From this we extract the dynamical density of states at 
finite temperature and look for sign of a gap. 

Introduction 



Motivated by a series of recent experiments[l], [2], [3], [4], [5] Tsai et.al[6] have 
studied the question of reproducing the results of Eliashberg's strong couphng 
theory of superconductivity without making an assumption of a broken symme- 
try phase. The need for this stems from fact that in real systems, there may be 
many sources of instability and there is no general guide (other than through 
experiments) to determine which of these instabilities dominate. Thus from a 
theoretical standpoint, as Tsai et.al. [6] argue, it is desirable to have a theory 
that is unbiased in the sense that it makes no assumption about the system 
being in a given phase. Rather all sources instability are treated on an equal 
footing and a sufficiently powerful theory should be able to pick out which of 
these phases dominate in what regions of the coupling constant/temperature 
plane. The RG approach of Shankar[7] as applied to the problem of strongly 
coupled superconductors by Tsai et.al. [6] is one such. Here we present an 
alternative to this interesting and important work, namely the hydrodynamic 
approach. In the hydrodynamic approach, we simply recast the electron-phonon 



+ electron-electron problem in terms of the electron's hydrodynamic variables 
namely the density fluctuations and the conjugate namely the velocity potential. 
The crucial fermionic nature of the electron is captured by a nontrivial phase 
functional that enters into the framework in the lagrangian formulation. The 
hamiltonian is identical to the Dashen-Sharp formula in terms of currents and 
densites. 

Unlike the RG approach^, which requires dividing the Fermi surface up into 
patches, we do not have to pay conscious attention to the Fermi surface. The 
properties of the free Fermi theory are automatically encoded in the phase func- 
tional. However our approach is by no means exact. The phase functional that 
encodes Fermi statistics has to be determined by expanding in powers of den- 
sity fluctuations. We retain only the leading linear term, thereby implying that 
three-body density correlations are implicitly ignored. This means we have to 
work in the high density limit, where the Migdal's theorem applies. This is 
exactly the regime in which Tsai et. al. operate. However, our approach can be 
generalised in principle to regimes of lower density as well by including three- 
body terms and so on (but the calculations are obviously harder to carry out 
in practice). Thus the purpose of the present article is to highlight the useful- 
ness of the hydrodynamic approach by comparing favorably with the results of 
well-established methods authored by famous physicists^. 

2 The Hydrodynamic Action 

In an earlier preprint [8], we recall that we had written down the action for 
electrons coupled to phonons in the hydrodynamic language. We rederive that 
here for the sake of completeness. We assume following Tsai et.al.[6] that the 
electron-electron interaction is short-ranged as is the electron-phonon interac- 
tion. In their work, the bare electron-electron coupling constant is denoted by 
«o and the bare electron-phonon coupling g. We try to follow the notation of 
Tsai et.al. as closely possible. The full action may be written down as follows. 
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The prime over the summation q indicates that we restrict the values of |q| < Ad 
the Debye cutoff. The free part may be recast in the hydrodynamic language 
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as follows. 
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Here the phase function (f>e([/9];x(T) encodes Fermi statistics. Without it, we 
would be describing bosons rather than fermions. Also the Fermi current is given 
by the hydrodynaniic expression which implicitly defines 11 namely, J = — pVH. 
The phase functional $e has to be fixed by making contact with the correlation 
functions of the free Fermi theory. This has been done in an earlier preprint and 
we may write down the leading term in the expansion in powers of the density 
fluctuations. 

$e(M;xa) = ^e'l-e^"* C{q,na) p_q,,_„ (3) 

qn 

Here N'^ is the total number of electrons and, 

[Pcia{t)p-cia{t))^ = Y, e^"'*^* ^ {Pqa.nP-qa.-n)o (5) 
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is the density correlation function of the free Fermi theory. The slow part of the 
field variable is given by. 

Also Zn — 2^ is the bosonic Matsubara frequency. Furthermore, the myste- 
rious factor of two in the exponential is to ensure that the exponents in the 
one-dimenional case come out right. By expanding in powers of the density 
fluctuations and retaining the leading terms we may write the full action as 
follows. 
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Here pq_„ ~ Pq|,n + Pqi,n- FoUowing Tsai et.al. we choose the phonons to be 
Einstein-like namely constant dispersion Wq — we- Since the action is purely 



quadratic in the phonons, we may integrate it out and write down an effective 
action for the electrons. 
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Now we wish to compute the propagator. More specifically, we wish to simply 
compute the momentum distribution and from that extract the quasiparticle 
residue at finite temperature. Tsai et. al. use this procedure to ascertain the 
transition temperature for the metal-superconductor transition. However the 
crucial point is that the vanishing of the quasiparticle residue is a necessary 
but not sufficient condition for the system to be a superconductor. It could 
also be an insulator. Thus either one has to demonstrate the divergence of 
the d.c. conductivity or more illuminatingly demonstrate phase coherence. In 
other words the nonvanishing nature of Yang's offdiagonal long-range order 
correlation function. However since we know from prior experience, that this is 
a superconducting transition, we shall content ourselves in just computing the 
momentum distribution. After integrating out the phonons we may write, 
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The full propagator may be written as follows. 
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Now we make use of the trick outlined earlier namely multiply and divide by 
the free propagator, in the denominator use the the bosonized version and in 
the numerator use the one obtained from elementary considerations. Thus we 
may write (see appendices for definitions notation e.t.c), 
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(12) 
Thus we have an asymptotically exact formula for the propagator provided we 
evaluate E. To this end we first integrate out the down spin variables and write, 
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We may now read off the correlation functions. 
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We use the current algebra constraint to determine the approximate long- 
wavelength nature of the correlation functions of the free theory. This means the 
unknown C has to be fixed so that the current-current correlation arc related 
in the usual manner to density-density correlation functions. 

JqT,« = ^-qT,-n(«— q) (22) 
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This in turn means, 
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Just to verify that this is sensible, we compute the static density-density corre- 
lation at zero temperature which we know is {N^ /2)q/{2kF)- 
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as required. It is important to point out that there is no such simple connection 
between density-density and curcnt-current correlation functions for interacting 
systems, since the four-point functions are not obligated to resemble the non- 
interacting values that enable the correspondence. Therefore, in particular one 
should not look to similarly relate Eq.( 16) and Eq.( 17). Finally, we wish to 
ascertain that the last correlation function is consistent with current algebra. 
We know that. 
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Thus we must have, 

an identity. Now we may proceed to evaluate the fuU propagator. 



3 Dynamical Density of States 

Here we compute the full propagator (see appendices for more details and def- 
initions of the various terms). From there we extract the dynamical density of 
states. An examination of this should tell us whether or not a gap is persent. 
The one-particle spectral function is given by, 

2TT A{Ku;) = y" dt e'-* [{c^^{t)cl^{0)) + (4^(0)ckT(t))) (32) 

We define the dynamical density of states as, 

k 

Thus the relevant quantities are the unequal-time, equal-space Green functions 
: ( T ip{x |,t)'0t(x "fjt )). From the appendices we find that we may simplify 
these in two spatial dimensions. 
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To first order, we take the zero temperature free Green functions. 
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as required. In general, 
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where, 

A = fc^T ^(V/^^ - ^o)^D (43) 
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We may provisionally identify A with a gap. Thus we find that there is no 
superconductivity unless the (screened) phonon strength exceeds the (screened) 
electron-electron repulsion. Thus wc have the (necessary) condition for super- 
conductivity g^ > uoWe/2. In particular we may also conclude that for g^ = 
and Mo < also we have a gap and hence possibly also superconductivity. 
Strictly speaking we have to compute the d.c. conductivity also and show that 
it diverges in order to be convinced that the transition is superconducting (or 
demonstrate phase coherence ). But since we already know this to be a super- 
conducting transition we shall content ourselves with doing the bare minimum 
as do Tsai et.al. 

This preprint is Incomplete 

I seem to be making some errors in the computation of 
integrals. Perhaps some more knowledgeable people can 



help me out. In any case this is meant to show members 
of hiring commitees that I am despertaely trying to do 
something important. 

4 Appendix A 

Here we provide some details of the computations of the propagator. 
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We note that the pathological n = should be excluded from consideration. The 
reasons for this are not entirely clear but doing so enables the right exponent of 
the Hubbard model to be recovered[9]. When m — 1,2 we have, 
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In the noninteracting limit, this sum is zero since w* -^ vp and w^ -^ kp and 
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5 Appendix B 

Here we calculate the equal-space unequal-time Green functions. 
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The niomentum integrals cannot be done exactly, hence we make some assump- 
tions. We assume that the critical temperature lies betwen the Debye energy 
and the phonon energy. In other words, vpAjj <C ksTc <C we- This means, 
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First we focus on two dimensions. We have to make use of the temperature 
constraint repaeatedly namely vpAjj 't^ ksT '^ we- Then, 
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Since t — t ~ l/lksT) the temperature constraint tells us that we may expand 
in powers of small q. 
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